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^ ■ Abstract 

We provide exact analytical expressions for the magnetic susceptibility func- 
tion in the high temperature expansion for finite Heisenberg spin systems 

^^ ' with an arbitrary coupling matrix, arbitrary single-spin quantum number, 

Q ' and arbitrary number of spins. The results can be used to determine un- 

known exchange parameters from zero-field magnetic susceptibility measure- 
ments without diagonalizing the system Hamiltonian. We demonstrate the 
possibility of reconstructing the exchange parameters from simulated data for 
two specific model systems. We examine the accuracy and stability of the 
proposed method. 
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I. INTRODUCTION AND SUMMARY 



m 



Due to advances of synthesis methods in coordination chemistrya'B and, especially, 
polyoxometalate chemistry,ErQ the field of molecular magnets has undergone rapid progress 

Ch ! in recent years. To date each of the identical molecular magnets contains as few as two and 

up to several dozens of paramagnetic ions ("spins"). The largest paramagnetic molecule 
synthesized, the polyoxometalateQ {Mo72Fe3o} contains 30 iron ions of spin s = 5/2. Al- 
though these materials are synthesized as macroscopic samples, i.e. crystals or powders, 
the intermolecular magnetic interactions are generally utterly negligible as compared to the 

5t ! intramolecular interactions. Therefore, measurements of their magnetic properties reflect 

mainly ensemble properties of single molecules. It appears that in most cases the magnetic 
properties of these molecules are rather well described by the Heisenberg model with possibly 
small anisotropy corrections.lIrEj 

In order to calculate from first principles the energy spectrum and the magnetic prop- 
erties it is essential to accurately determine the exchange parameters, as these measure the 
strength of the two-spin interactions. The usual procedure consists of the following steps: 
diagonalization of the system Hamiltonian for an arbitrary choice of exchange parameters; 



calculation of the partition function and a quantity, such as the magnetization, which can 
be compared to experimental data; variation of the exchange parameters so as to achieve an 
optimal fit between the experimental and theoretical data. This can be a rather straightfor- 
ward task for simple configurations such as small spin rings with a single nearest-neighbor 
coupling. However, for more complicated systems, for example where two or more exchange 
parameters may be needed it may be impractical to perform such calculations for different 
trial choices of the parameters, even if one is guided by chemical and structural information. 
Worse yet, for large spin systems, where in practice the Hamiltonian cannot be diagonalized, 
it may even be impossible to fix the exchange parameters, since one cannot calculate the 
weak-field susceptibility over the full tempearature range, for comparison with experiment. 

In this article we explore the possibility of determining the exchange parameters of 
Heisenberg molecular magnets using the exactly calculable first three terms of the high- 
temperature exDansion of the zero-field susceptibilty xiP): iii powers of (3. In comparison 
to earlier workJl^EJ exact results are derived here for these three coefficients, for arbitrary 
coupling, arbitrary single-spin quantum number, and an arbitrary number of spins. The 
procedure we explore here opens up the possibility of treating large systems where the 
Hamiltonian can no longer be diagonalized. We also outline the necessary algorithms for 
calculating additional terms of the high-temperature expansion, but these would require 
computer algebra methods. 

Deducing exchange parameters from experimental susceptibility data of course presumes 
that one possesses accurate estimates of the number of paramagnetic ions, the spectro- 
scopic splitting factor (s), diamagnetic contributions, temperature- independent paramag- 
netism, etc. We do not consider these issues here; instead we concentrate on the question, 
with what accuracy exchange parameters can be deduced using theoretical susceptibility 
data, for a Heisenberg molecular magnet, limited to the high-temperature regime. We do 
not want to mix questions about the accuracy of the method with uncertainties in experi- 
mental data which are of different origin. 

This article is organized as follows. After recalling basic definitions in Sec. |I| we derive 
the trace formulae which are necessary for the high-temperature expansion in Sec. |T|. These 
results are summarized in Sec. |IV| and related to the expansion coefficients of the Taylor 
series of xiP)- The accuracy and stability of the proposed method are illustrated for two 
selected model systems in Sec. |V|. We briefly summarize our flndings in Sec. |VI[ 

II. DEFINITIONS 

For the following investigations we assume that the spin system under consideration is 
well described by the following general bilinear Heisenberg Hamilton operator 

3 

¥o = E ^m^Im ■ i^ = E ^m^ E ^si^^ ■ (1) 

The exchange parameters J^^, are given in units of energy, with J^^ < resulting in ferro- 
magnetic, J^i, > in antiferromagnetic coupling, s^ are the dimensionless single-particle 

spin operators and Sjl denotes the spin components i G {1,2,3} = {x,y,z} of the spin 



at site /i. The interaction with the magnetic field is taken into account by including the 
Zeeman term 

H = Ho + 9I2bBS^'^ , (2) 

where S"*-^-* is the 2;-component of the total spin operator, and g is the spectroscopic splitting 
factor. For simplicity we are supposing that all of the paramagnetic ions possess the same 
value of g. The present results can be easily generalized to cases where several different 
paramagnetic species are present. As usual the magnetisation is defined as 

-M = -^Tr|^/iB§(^)e-^^} , Z = Tr {exp(-/?^)} , (3) 

Z being the partition function. Accordingly the zero-field magnetic susceptibility is 

fdM\ Tr{g(3)2exp(-/5^o)} 
Xo = ^ = (3{g^^Br^-^- ^ . (4) 

The susceptibility may be renormalised by division of [gfiB)^ 

A Taylor expansion of x in powers of /?, starting from (^), yields after some calculations 

m = f: c„r = ^P- ^P' + {--^ + ^) P' (6) 

n=l ^0 ^0 V 



'^\6ty 2tl 6tJ^ 

+ ^^°U^ 24tlJ 6tl 4t2+24tc'^ 

+ ... , 



where 



Comparing expansion @ with a Curie- Weiss law 

XcwiP) = Y&f = ^^^^^^ ' ^^'^^'^ = ^^^^ ^^ + ^^^^^ + {kBT,l3f + ...) (8) 

yields the known theoretical resultsij for the Curie constant C and the Weiss temperature 
% 

ksC = ^ and kBT, = -^ . (9) 

to /io 

Generally, the /3^-term of (|^) will not coincide with the /5^-term of the rigorous expansion (^ 
This is due to the fact that the Curie- Weiss law itself is only valid at high temperatures.li3 
Thus, the /5^-term of (^ constitutes a first (quantum) correction to the Curie- Weiss law for 
lower temperatures. 



III. TRACE FORMULAE 

Expanding the terms occuring in the trace (|^ 

t„ = Tr jff^} = Tr I 5^ J,. ^ .«. W 1 (10) 

and using Tr(y4 ® B) = (Try4)(Tri?) one ends up with terms of the form 

Tr(A...A,),A,G{s«,£(2),,(3)}. (11) 

Here the spin operators s without a site index denote operators operating in the single-spin 
Hilbert space C^*^^. Let (ii{i = 1,2,3) denote the number of occurrences of s*-*-* in the 
product Ai . . . Ai. One can easily show that Tr(y4i . . . Ai) is non-zero only if all li are even 
or all li are odd. We give a list of the simplest cases, where the trace is non-zero: 

£ = 0:Tr{l^} = 2s + l, (12) 

£ = 2 : Tr js^^^'} = ^ m^ = -s{s + l)(2s + 1), 

m=—s 

£ = 3 : Tr {^^'^s^'^s^'^} = -Tr {s^'^s^'^s^^^} = '-s{s + l){2s + 1), 

s , 

£ = 4 : Tr js^^} = ^ m^ = —s{s + l)(2s + l){'is^ + 3s - 1), 
£ = 4 : Tr jsW's^')'} = —s{s + l)(2s + l)(2s' + 2s - 1) . 



Note that there are no £ = 1 cases; hence all traces Tr < s}j\ . . . Sj^k \ vanish if at least 
for some k = 1, . . . ,K the site index /i^ is different from the other //'^ and jk = 1, therefore, 

ti = Tr|^o}=0. (13) 

From the list (|T2D we obtain for n = 2: 

t, = Tr {hi} =Y: Jl E Tr {s^^ ® s^' ® 1^_,} (14) 

where 1 n-2 denotes the unit operator in the Hilbert space of A^ — 2 spins. For n = 3 there 
occur two kinds of non-zero terms: 

-^{s^s^s^^^s^'s^^In^^ (15) 

and those terms obtained by permutations of {1, 2, 3}, and 



Trj^f ®£«^®£«^®liV-3) , (16) 
where /i, u, k are pairwise distinct and i = 1,2, 3. Consequently we obtain 

t3 = Tr|^3} =Ti + r2, (17) 

Ti = ^ J3^ 3!(2s + l)^-2Tr {s^'^s^'^s^^^^}' (18) 

T2 = 3! [ 5^ J,. J,. J. A (2s + 1)^-3 X^Tr{£«2|^ (^9) 

\fJ,<l^<K / 



IV. THE TAYLOR EXPANSION OF x 

A. The linear term 

For the hnear term we just need to and fiQ. to is nothing else than the dimension of the 
A^-spin Hilbert space 

^o = TrU} = (2. + l)^. (20) 

For fiQ we obtain with the help of list ([l2D 

Tr|s(=^)2|^i.j.^||2| (21) 

I (ETr{l'}+ETr{|.-|M}) =lNs{s + l)to 



yWo 



3 



f-O 



and therefore 



ci = ^Nsis + 1) . (22) 



Hence, combining (^, (pOD, and (^) we find that the Curie- Weiss constant is given by 

knC = kBigfisfC = [gfi^y^Nsis + 1) , (23) 



which is identical to the mean-field result .HjEJ 



B. The /32-term 

For the quadratic term we calculate /xi directly 

^,, = Tt{HoS'^'^'] = ^-Tt{HoS'} (24) 

= I E ^M^Tr L«.« (2 Y: isi^^s!^^) + Ns{s + 1)^] I . 

IM<U,i (. \ K<X,j / J 

Since Tr{sfi sj } = for /i 7^ z/ the second summand vanishes. According to section |T|, 
the first one is only non-zero ii n = k,!/ = X,i = j. Hence 

/^i = ^ E ^/^^Tr {.f £«^} = 2 [5^ J,.] Q.(. + 1))' (2s + 1)^ (25) 

fj.<u,i Vm<'^ / 
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This results in a coefficient 



C2 



-^(E-^-) -'(- + !)' (26) 



and a Weiss temperature of 

/ii 2 



which agrees with the mean field result .Ej'Ej 



kaTc = = s(s 

/io 3A^ 



C. The /3^-terni and higher orders 

Using the same method the next term can be calculated. The coefficient reads 
2 

C3 = —S%S + lf J2 {JnxJKt. + JKt.Jx^ + J.xJxt.) (28) 

K<A</1 



y{s + ifY.iJl>) ■ 



18 

K<A 

Hence C3 depends not on all exchange parameters but only on a certain combination of the 
second moments (or correlation coefficients) of the exchange parameters. 

It is clear how coefficients of higher orders are to be calculated, although the expressions 
rapidly become very lengthy. Computer algebra methods would have to be employed to 
determine and evaluate those expressions. 



V. APPLICATIONS 

The first two expansion coefficients in the Taylor expansion (^ of the susceptibihty in 
powers of P that involve the exchange parameters are C2 and C3. We propose exploiting ac- 
curate numerical estimates for C2 and C3 to provide information on the exchange parameters. 
The following version of (P) 

00 
kBT{kBTx-ci) = J2cnf3''-^ = C2 + C3(3 + --- , (29) 

n=2 

is especially convenient since we can determine C2 and C3 by a least-squares fit between a 
low degree polynomial in /? to measured values of the quantity ks T {ks Tx — Ci), compare 
also Ref. 0. To increase the accuracy of the determined values of C2 and C3 it is generally 
advisable to use a polynomial in (3 of degree two or higher. In an iterative procedure one 
could check the extent to which the values of C2 and C3 change as higher degree polynomials 
are used. We remind the reader that we assume that the spectroscopic splitting factor g, 
which is necessary to obtain function (|29|) , has accurately been determined, e.g. by ESR 
measurements. 

In the following two examples we test the proposed method for two model systems. 
Starting with given exchange parameters we determine the partition function by exact di- 
agonalization of the Hamiltonian and calculate the susceptibility. Our aim is to extract 
accurate estimates for the exchange parameters from a limited number of susceptibility data 
points. We also assess how the method depends on the chosen temperature interval and how 
robust the method is when the experimental data includes random error. 

If the system has two exchange parameters Ji and J2 it is obvious from (p6| ) and ( |28D 
that C2 is a linear function of Ji and J2 whereas C3 is a quadratic function of these variables. 
Hence Ji and J2 can be determined by numerically solving a simple non-linear equation once 
accurate values of C2 and C3 are available. In cases where multiple solutions exist additional 
magneto-chemical information, for instance concerning the sign of the exchange parameters, 
must be considered. 

A similar complication arises for the special class of so-called isospectral spin systems, 
where the same eigenvalue spectrum is associated with a host of different choices of the ex- 
change parameters.ll3 For such systems it is in principle impossible to obtain unique exchange 
parameters from susceptibility data. 

A. Cubic spin array 

Our ffist model system is a spin array where A^ = 8 spins s = 1 are located at the 
vertices of a cube. The spins are assumed to interact via two antiferromagnetic exchange 
parameters of typical magnitude, Ji/kB = 20 K between nearest neighbors in the top and 
bottom face and J2/^s = 10 K connecting the two faces, see Fig. |1|. The calculated magnetic 
susceptibility is shown in Fig. ^ as a thick line. The thin horizontal line corresponds to the 
high-temperature limiting value of ksTx given by the coefficient Ci, see (|^) and (P^. For 
the determination of the underlying exchange parameters a restricted set of susceptibility 



data is used. This set consists of data taken in steps of 5 K from 150 K to 300 K (crosses 
in Fig. I). As already mentioned, not the susceptibihty x itself or Tx are fitted but 

f{T) = T{kBTx-c^). (30) 

The graph of the latter function is shown in Fig. 0. One notes that its dependence on 1/T 
in the relevant temperature interval is already rather linear. Therefore one expects good 
results using the fitting procedure based on a low degree polynomial. Indeed, approximating 
T [ksTx — Ci) by a polynomial of third order yields Ji/ks = 19.96 K and J2/kB = 10.09 K, 
a result which is not improved by using polynomials of higher degree. 

We also investigated how the chosen temperature interval influences the analysis. In- 
cluding only data between 150 K and 200 K yields exchange parameters which deviate by 
at most 10 % from the correct values. 

Different from our model investigation, experimental data necessarily incorporates ran- 
dom noise. In order to estimate this effect, the determination of the exchange parameters 
was carried out with a number of modified data sets. For each data set all susceptibility 
values were allowed to deviate by at most ±5 % from the correct values. The resulting 
exchange parameters are shown in Fig. ^, each point corresponding to a different modified 
data set. It appears that the influence of the fluctuations is different for the two exchange 
parameters due to the non-linear dependence of Ji and J2 on the expansion coefficients c„. 

B. Ring spin array 

As a second model system we investigate a ring with A^ = 6 and s = 5/2. In order to 
further quantify the accuracy of the proposed method, we chose two cases, case (a) where 
the spins are connected only by a nearest-neighbor interaction Ji/ks = 20 K, and case (b) 
where we have an additional next-nearest-neighbor interaction J2/^_b = 10 K, see Fig. ^. The 
corresponding values of Tx, again obtained by complete diagonalization of the Hamiltonian, 
are given in Fig. |. The thick curve corresponds to case (a) and the dashed curve corresponds 
to case (b). The thin line again gives the high temperature limit. Data points used for our 
analysis are depicted by crosses, again within the range between 150 K and 300 K. Note 
that for this model system the data points deviate from the high temperature behavior much 
more than in the previous example, i.e. the dependence of f{T) on 1/T (Fig. ^ in the above 
temperature range especially for case (b), is not as linear as it was in the previous example. 

The present method is sufficiently robust to provide accurate estimates of the exchange 
parameters without exploiting any prior information. This is illustrated by using two ex- 
change parameters for case (a) and finding that Ji/ks = 20.04 K and J2/^b = —0.04 K, 
which is in very good agreement with the actual exchange parameters in case (a). For case 
(b) we find Ji/ks = 19.33 K and J2/^_b = 10.69 K, which is not as accurate as in the 
previous examples. Nevertheless, even in this worst case, using these exchange parameters 
results in ksTx which, when plotted in Fig. ^, is indistinguishable for all temperatures from 
that for the correct values of Ji and J2. Improved exchange parameters might be achievable 
if measured data for x is available for higher temperatures, but unfortunately many sub- 
stances - especially those containing solvent molecules for instance - become structurally 
unstable at such temperatures. 



VI. SUMMARY 

In summary, we have provided exact analytical expressions for the first three terms of 
the high-temperature expansion of the magnetic susceptibility function for Heisenberg spin 
systems with an arbitrary coupling matrix J^j, and arbitrary intrinsic spins s. The results 
can be rather easily used in order to determine unknown exchange parameters for molecular 
magnets from zero-field magnetic susceptibility data with good accuracy, thereby obviating 
the need for diagonalizing the system Hamiltonian. 
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FIG. 1. Structure of a spin cube of N 
Ji/kB = 20 K, J2//CB = 10 K. 



and s = 1. The assumed exchange parameters are 
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FIG. 2. ksTx data (thick hne) for the spin cube shown in Fig. |l[ The data set used in the 
calculations is denoted by crosses. The thin line is the high temperature limit. 
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FIG. 3. Functional dependence of T {ksTx — ci) on 1/T (thick line) for the data set (crosses) 
shown in Fig. |2[ The continuation of this function for higher temperatures is given by the dashed 
line. 



11 



12 



10 



>^_ 



19.0 



195 



200 20.5 

]l (K) 



21.0 



FIG. 4. Distribution of exchange parameters of the spin cube under the assumption that the 
susceptibihty data incorporates random errors at the level of 5 %. 




FIG. 5. Structure of a spin ring of A^ = 6 and s = 5/2 with nearest-neighbor Ji and 
next-nearest-neighbor coupling J2. The assumed exchange parameters are (a) Ji/ks = 20 K, 
J2 = and (b) Ji/ks = 20 K, Ja/A;^ = 10 K. 
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FIG. 6. ksTx data for the spin ring introduced in Fig. p. The thick curve shows case (a) and 
the dashed curve shows case (b). The thin line corresponds to the high temperature limit. The 
data sets used in the calculations are depicted by crosses. 
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FIG. 7. Functional dependence of T (/c^Tx — ci) on 1/T (thick line) for the data set (crosses) 
shown in Fig. |y. The continuation of these functions for higher temperatures is given by dashed 
lines. 
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